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We argue based on theoretical considerations and analysis of experimental data that quasiparticle excitations 
near the nodes determine the low temperature properties in the superconducting state of cuprates. Quantum 
effects of phase fluctuations are shown to be quantitatively important, but thermal effects are small for T <^Tc. 
An anisotropic superfluid Fermi liquid phenomenology is presented for the effect of quasiparticle interactions on 
the temperature and doping dependence of the low T penetration depth. 



1. INTRODUCTION: 

The superconducting (SC) state of the high Tc 
cuprates differs from conventional SCs in several 
ways: a d-wave gap with low energy quasiparticle 
excitations near the nodes, a small phase stiffness 
and a short coherence length. There is some con- 
troversy about the importance of quasiparticles 
versus phase fluctuations j2j in determining 
the low temperature properties. In this paper, we 
discuss this problem focusing mainly on the dop- 
ing and temperature dependence of the in-plane 
superfluid stiffness D^^ which is related to the pen- 
etration depth A|| through = A-Ke^D^Jfi^c^dc 
where dc is the mean interlayer spacing; we will 
set Ti = c = e = \ below. 

We first review experimental evidence for 
quasiparticle excitations at optimal doping. 
Transport data in the SC state in YBCO [||,|) 
shows a scattering rate decreasing sharply be- 
low Tc implying long lived quasiparticle excita- 
tions for T <t: Tc. Direct evidence from ARPES 
in Bi2212 shows the presence of sharp quasipar- 
ticle peaks over the entire Fermi surface for 
T < Tc. Thermal conductivity data §] in YBCO 
and Bi2212 shows k ~ T at low T. The slope 
predicted by quasiparticle theory |Q is in good 
agreement with this k data on Bi2212, using 
ARPES Q estimates for the Fermi velocity vp 
and the gap slope wa — {2hkp)~^dA/d(l) at the 
node. Experimentally, it thus seems that quasi- 
particle excitations exist and are important at low 



temperature. 

It then seems natural to interpret the linear 
T dependence of A|| (T) Q as arising from nodal 
quasiparticles (QP). Ignoring QP interactions the 
layer stiffness D^^ (T) = (0) - AqT with Aq = 
{kgln2/TT)vF/vA. ARPES estimates in Bi2212 
for vf and wa give ^ 0.8meV/K, whereas 
experiments at optimality measure a slope ^ 
0.3 — O.AmeV/ K . Thus there is at least a factor 
of two discrepancy which needs to be understood. 

Alternatively, it has been suggested that this 
linear T behavior could arise entirely from ther- 
mal phase fluctuations |10 1^ without invoking 
nodal quasiparticles. However, there are two rea- 
sons to believe that thermal phase fluctuations 
are unimportant at low T in the cuprates. (1) An 
effective action calculation for charged d— wave 
SCs 1 12 1, summarized below, shows that thermal 
phase fluctuations become important only near 
Tc for Bi2212 at optimality. On the other hand, 
quantum phase fluctuations are important at low 
T and suppress both D^^ (0) and the slope. For 
Bi2212, the resulting renormalized D^^ is about 
30% smaller while the renormalized slope is about 
25% smaller than the bare values. (2) Further, 
with underdoping, D^^ (0) decreases Q and the 
slope of Z?|| (r) also shows evidence of decreasing 
in Bi2212 and La214, although some YBCO data 
is consistent with a doping-independent slope 
(see the compilation in refs. ^,|l^). Insofar as 
the data indicate a doping dependent slope for 
-D|i (r), they independently rule out classical ther- 
mal phase fluctuations as the explanation pO| , pl| 
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for the linear T dependence, since the slope of 
Z?ii in such theories is insensitive to doping. Fi- 
nally there is another proposal for understanding 
{x; T) invoking incoherent pair excitations ||l5[| 
which, in our opinion, does not properly include 
the effect of Coulomb interactions. 

We next summarize our phase action calcula- 
tion, and then describe how quasiparticle inter- 
actions could account for the difference between 
the free QP value and the measured slope and its 
doping dependence. 

2. PHASE FLUCTUATIONS: 

We have recently derived, by appropriate 
coarse-graining, a quantum XY model describ- 
ing phase fluctuations in charged, layered d— wave 
SCs, starting with a lattice model of fermions. 
See ref. [|l2| for details of this derivation and some 
of the discussion in this Section. 

The phase action for layered SCs with in-plane 
lattice spacing a = 1 takes the form 

1 /■^/'^ F 

+ - / dr^D^ [l-COs(6lr,r-6'r+a,r)] 

where is the in-plane coherence length, 

and D refers to the layer stiffness without 

II ■' 

phase fluctuation effects, but including possi- 
ble renormalizations due to quasiparticle inter- 
actions. Here — V"(q|| /fo, qj_) with V{c\) = 
(27re2/q'||e6)sinh((7||4)/ [cosh(g||(ic) - cos(q_L4)] 
is the Coulomb interaction for layered systems, e?, 
is the background dielectric constant, dc the in- 
terlayer spacing, and Qh , q^^ refer to in-plane and 
c-axis momentum components. The prime on the 
sum indicates a Matsubara frequency cutoff since 
the energy of the fluctuations should not exceed 
the condensation energy i^cond = \D^^ (tt/^q)^- 
The form of the first term of (1), which arises 
from coarse- graining up to a scale ^Oi E^nd the 
importance of cutoffs have not been appreciated 
earlier. The (typically very small) c-axis stiff- 
ness D ^ can be ignored for in-plane properties 



since it was found not to lead to qualitative or 
quantitative changes. 

We ignore vortices (transverse phase fluctua- 
tions) which are suppressed at low T by their fi- 
nite core energy. Analyzing longitudinal phase 
fluctuations for (1) within a self-consistent har- 
monic approximation (SCHA) leads to the 
renormalized stiffness D^^ = D^^ exp{—{56'^)/2) 

where 66^ = {9^^- — 6'r+Q,r)^- Our numeri- 
cal results can be simply understood as follows; 
(<56'2)(r = 0) ~ ^(eVeb^o)/-D||(0) is a measure 

of zero point quantum fluctuations, while classi- 
cal thermal phase fluctuations become important 
near a crossover scale Tx ^ min [Tc,TO] where 

T° = y'i:'||(0)(e2/ebCo) is the T = oscillator 
level spacing in the renormalized harmonic the- 
ory. 

It is easy to see that phase fluctuation effects 
are negligible in the BCS limit, except very close 
to Tc. With eVefofl ^ d'^^ ^ and Co ~ ^f/A, 
one obtains the standard result E'cond ~ A^/i?^ 
per unit cell, and {59'^) [T = 0) - y^A/E^ < 1 
and Tx - min [T„ /E^A] = T^. 

F 

For the cuprates, the small and small D^^ 

act together to increase {S9'^), but they push Tx 
in opposite directions. For optimal Bi2212 we 
use e^/cha « 0.33ey with Ch ~ 10, <^o/a ~ 10) 
and dc/a ~ 4. Assuming that the two layers 
within a bilayer are phase-locked, we get the bi- 
layer stiffness D^^ (0) « %QmeV , from experimen- 
tal data which shows A,, (0) w 2000^. This leads 
to Scond ~ eif/unitcell and « QOOK. Since 

F 

the bare stiffness D^^ actually decreases with tem- 
perature due to quasiparticle excitations, an es- 
timate of the crossover scale Tx can be obtained 

from Tx ~ \l ^\\ (^x)(e^/ef,Co)- Assuming a lin- 
early decreasing D^^{T), this leads to Tx ~ Tc- 
Thus longitudinal thermal fluctuations are clearly 
unimportant at low temperatures. Quantum fluc- 
tuations are important at low temperatures since 
{59'^){T = 0) ^ 1 at optimality and detailed cal- 
culations |l2j lead to a 30% decrease of D^^ (0) 
and a 25% decrease in the slope. 

While it might appear that there could be a low 
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T crossover to thermal phase fluctuations due to 
the low energy c-axis plasmon (~ IK for Bi2212) 
in the anisotropic layered SCs, the phase space for 
these low lying fluctuations is too small to lead to 
a linear T behavior . Even in a purely 2D sys- 
tem with a low lying .ygjj" plasmon, the decrease 
in the phase stiffness due to phase fluctuations 
only goes as a large power law (~ T^). Quasipar- 
ticles are thus crucial in obtaining the observed 
linear temperature dependence. 

We next turn to the doping dependence of 
phase fluctuations. The (amplitude) coherence 
length ^0 is crucial in determining the effect of 
phase fluctuations. Since is determined by the 
pairing gap which appears to remain flnite as we 
underdope, we do not expect singular behavior 
in the phase fluctuations arising from the doping 
dependence of ^o- In this case, the dominant dop- 
ing dependence to phase fluctuations arises only 
from the singular behavior of the bare parameters 
in the phase action on underdoping. This sin- 
gular x-dependence in D^^ {x) is most naturally 
explained by quasiparticle interaction effects as 
discussed below. 

To estimate the doping dependence of phase 
fluctuations, we note that the core energy will 
lead to vortices being exponentially suppressed 
at low r, even as we underdope. Using the ex- 
perimental input Tc ~ Z?!! {x]T — Q) ^ X and 
assuming a doping independent longitudinal 
fluctuations within the SCHA lead to ^ 
and hence r° » T^. Thus, Tx ~ min [T^TO] - 
Tc, which implies that thermal phase fluctuations 
are unimportant at low T as one underdopes. 
Further, within the SCHA, {50^) - x-^'^, which 
would lead to a destruction of superconductivity 
at small enough x. However, we do not expect 
the SCHA to be valid close to this transition. 

3. QUASIPARTICLE INTERACTIONS: 

The increasing importance of interactions with 
underdoping is evident: D^^ {x\ 0) x [|l3| and 
the quasiparticle weight diminishes |l6| as one 
approaches the Mott insulator. We thus explore 
the possibility that residual interactions between 
the quasiparticles in the SC state can account for 
the value and doping dependence of the slope of 



D^^[x;T). To this end we use a phenomenologi- 
cal superfluid Fermi liquid theory (SFLT) 0,|l8|. 
All available experimental evidence on the ground 
state and low lying excitations suggests that the 
correlated SC state in the cuprates is adiabati- 
cally connected to a d-wave BCS state. We thus 
feel that SFLT may be a reasonable description 
of QP interactions in the SC state at low T , even 
though this formulation makes reference to a (hy- 
pothetical) T = Q normal Fermi liquid in which 
SC is induced by turning on a pairing interaction. 
(One could argue that approaching the supercon- 
ducting phase from the overdoped side at T = 0, 
one obtains a normal Fermi liquid to SC transi- 
tion.) 

The bilayer stiffness after including QP inter- 
action effects is given by, -D^(r) = /^^^'"(O) - 
ap2{kgThi2/T:)vp/v/^ where a^^lSp are Fermi 
liquid renormalizations. We will constrain the 
Landau QP interaction function by demanding 
Pp ~ X, consistent with experiments and then 
determine the doping trends in . 

To compute a^, [3^, using a standard Kubo for- 
mula in the quasiparticle basis, it is convenient 
to shift the origin of the Brillouin zone to the 
(tt, tt) point and describe the hole-like Fermi sur- 
face of Bi2212 in terms of an angle (f>. The Lan- 
dau /-function is denoted by f{4>, 4>'). We deflne 
(O)^ = J^''dcl,kAcl^)0{cj))/[2n\vFm. We get 

= l + 47r((«f^,(0K,((/.')/(</.,(/.')))00V(4.)* 
from the diamagnetic response of the free en- 
ergy d'^SF/dAl to an applied vector poten- 
tial. The current carried by nodal quasiparti- 
cles is then renormalized by the factor y^a^ = 
1 + {vFx{<i))f{4'n,(l>))4>/ ['^VFx{4'n)\, relative to its 
non-interacting value, where the nodes are at 
0„ = (2n - l)7r/4 with n = 1 . . . 4. 

We expand /(0, (j)') = i^m,m' [cos{m(j) + 

■m'(j)') + cos{m'(j)+m'(j)')] in a set of complete basis 
functions [|l9|, where m, m' = 0, ±1, ±2, . . . with 
square lattice symmetry imposing m + m' — Ap 
with p = 0, 1, 2, . . .. In an isotropic system only 
p = Q survives and and are </>- independent. 
However, as emphasized in ref. [ p8| one then ob- 
tains = ~ in disagreement with experi- 
ments |9|. 

To illustrate how anisotropy can qualitatively 
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Figure 1. Doping dependence of renormalizations 
and pp plotted for anisotropic model discussed in the 
text. For the isotropic case one gets — Pp. 

change this scahng we keep only the leading p = 
term: f{(j),(j)') = 2Fi^i cos{(j) — <7i')i but retain the 
full anisotropy of the dispersion seen in ARPES 
poj . We make a reasonable choice of = P + 
Qx, with P such that (3^, ^ x as x ^ 0, and Q 
such that Pp — 0.5 at x = 0.2. This leads to 
ap{x) shown in Fig.l, which is a weak function 
of doping. In general there are too many free 
parameters in the anisotropic case (an infinite set 
Pm,m') for the theory to have predictive power; 
nevertheless the simple example above shows how 
the T = value and slope of D^^ can easily exhibit 
rather different x-dependences, and account for 
the experimentally observed D^^ (x; T). 

We thus arrive at the following picture for the 
doping and temperature dependence of D^^. The 
bare stiffness arising from non-interacting quasi- 
particles is renormalized by both QP interactions 
and quantum phase fluctuations at low T leading 
to the measured stiffness, D^^{x;T). Its doping 
dependence, D^^ {x; 0) ~ x, is determined by QP 
interactions while its linear T behavior is gov- 
erned by nodal QPs, with its slope renormalized 
by both QP interactions and quantum phase fluc- 
tuations. 
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